A New Method of the High Temperature Series Expansion 
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We formulate a new method of performing high-temperature series expansions for the spin-half 
Heisenberg model or, more generally, for SU(n) Heisenberg model with arbitrary n. The new 
method is a novel extension of the well-established finite cluster method. Our method emphasizes 
hidden combinatorial aspects of the high-temperature series expansion, and solves the long-standing 
problem of how to efficiently calculate correlation functions of operators acting at widely separated 
sites. Series coefficients are expressed in terms of cumulants, which are shown to have the property 
that all deviations from the lowest-order nonzero cumulant can be expressed in terms of a particular 
kind of moment expansion. These "quasi-moments" can be written in terms of corresponding 
"quasi-cumulants" , which enable us to calculate higher-order terms in the high-temperature series 
expansion. We also present a new technique for obtaining the low-order contributions to specific 
heat from finite clusters. 

keywords: high-temperature expansion, Heisenberg model, correlation function, specific heat, 
new method 
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I. INTRODUCTION 

The high-temperature expansion can be used in any 
dimension and has provided significant information on 
variety of models.EJ It is based on the Taylor expansion 
of the Boltzmann factor e~'^^ in /3, around the high- 
temperature limit. Although the concept of the high- 
temperature expansion is quite simple, its high-order cal- 
culation requires sophisticated strategies. A standard 
method of the high-temaerature expansion is called the 
finite cluster methodldo'Ll (FCM). It has been used for a 
long timeP and progress of the high-temperature series 
has been due mainly to improvement of computers. In 
the FCM, the series expansion is reduced to calculation 
in connected finite-size clusters. However, the FCM is 
inefficient for calculating a correlation function between 
rather distant sites: (i) The clusters have to include those 
two sites and calculation in such large clusters consumes 
a lot of time, (ii) The FCM does not use a valuable piece 
of information. That is, the high-temperature series of 
a long-range correlation starts at a certain finite order, 
namely, low-order coefficients are equal to zero phys- 
ically, it means that a long-range correlation develops at 
low temperature. 

In this paper, we formulate a new method oriented 
to the long-range correlation, namely, it overcome (i) 
and (ii) above. The new method is valid for the spin- 
1/2 Heisenberg model, or, more generally, for SU(n) 
Heisenberg model with arbitrary n. High-temperature 
series coefficients are written in terms of cumulants. We 
start from the lowest-order nonzero cumulant in the 



correlation-function series. Next, we consider a deviation 
from it. Then, the deviation can be regarded as a sort of 
moment, which we call a quasi- moment. The correspond- 
ing quasi-cumulant enables us to calculate a number of 
terms higher than the lowest-order nonzero cumulant. In 
Ref. ^, the Fourier transform of the correlation function 
is calculated using the new method and the FCM comple- 
mentarily. That is, we have used the new method for the 
long-range correlations and the FCM for the short-range 
correlations. 

In addition, also for the specific heat, we have adopted 
a similar strategy in Ref. ^. That is, contributions from 
the required largest cluster are calculated by a new tech- 
nique, and the FCM is used only for the smaller clusters. 
We utilize cumulants relevant to the ordering of quantum 
variables. By choosing nonzero cumulants from them, 
some of the series coefficients are calculated simply. We 
explain also this new technique in this paper. 

The structure of this paper is as follows. In Sec. 0, 
we present the model. We show a representation of 
SU(ri) Heisenberg model in terms of exchange operators. 
The relation of a correlation function and an expecta- 
tion value of the exchange operator is mentioned. In 
Sec. [II, we review the high-temperature expansion. Al- 



though we use this model for illustration, what is written 
here can generally be applied to other models. To clarify 
the purpose of our new method, we review the FCM in 
detail using a one-dimensional system as an explicit ex- 
ample. In Sec. IV, we describe the new method for the 
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correlation function. Here, we derive some properties of 
the cumulants. Owing to those properties, we can de- 
fine quasi- moments/cumulants. How to calculate series 
coefficients using those quasi-cumulants is explained. In 
Sec. 0, we show the new technique for the specific heat. 
Finally, we give a summary in Sec. VI. 
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II. MODEL 

Let us consider the spin-1/2 Heisenberg model, 

'^SU(2) ^ 5Z '^^J ' ^"'^■^ 

where Si is the spin operator at site i, and Jij is the 
couphng constant. Here, represents that the sum- 

mation is performed only once for each pair. The 
spin operators are rewritten as 

2 Si ■ Sj — Pi J — — , (2) 

where Pij is an exchange operator defined by 

••■ A... ^...) = I •■•^ •■•«••■), (3) 

for arbitrary a =1,1, and (3 =t, i. Then, we can rewrite 
Eq. §) as 

- E J^^^ (^^. - 1) ■ (4) 

This Hamiltonian has SU(2) symmetry. Since our 
method is valid more generally for an SU(n) symmet- 
ric case with arbitrary positive integer n, we define the 
SU(n) Heisenberg model below. 

Let each site take one of the n colors, and denote them 
as la) with a = 1,2, ■ ■ ■ ,n. An exchange operator is 
defined using := |q;)(/3| by 

n n 
a=l (3=1 

for i ^ j- Colors at sites i and j are exchanged by P^j, 
namely, Eq. is satisfied for arbitrary colors a and [3. 
Furthermore, for a later convenience, let us define 

Pi,] ■■= Pi,j - ^- (6) 
In the case of SU(2), 

=2 Si - Sj (7) 

is satisfied. Although we mostly use expressions for gen- 
eral n in this paper, we expect that Eq. helps instant 
understanding of explanations. In some cases, we explic- 
itly use the SU(2) notation for a simple explanation. 

The Hamiltonian of the SU(n) Heisenberg model is 
given by 

^^"^ ■■= E J^^^ (p^^^ - ^ E J^^^p^^^- (8) 

The constant term just yields an energy shift. Hence, for 
calculating thermal averages, we can use 

H:=Y,A,P^,r (9) 



We define the partition function as 

Z := Tr e^*^^, (10) 
and the free energy as 

F:=-p-HnZ. (11) 
The average in this system is denoted by 

(O) := Tr (O g-'^^) /Z, (12) 

where O is an arbitrary operator. 

In the case of spin-1/2 system, there are three indepen- 
dent interacting components, s^, s^, s^. In general, when 
there are n states per site, the maximum number of in- 
dependent interacting components is — 1. The model, 
isotropic with respect to these components, is the SU(n) 
Heisenberg model as explicitly shown in Appendix 

What we present in this paper is the high temperature 
expansion of {Pi,j)- A correlation function {Xf^X^"') is 
calculated using a relation, 

{Xfx^-) = ^ f (P,,> - i) , (13) 

for arbitrary a ^ (5, i ^ j . In addition, when i = j, 

{Xfxf'^) = i, (14) 
n 

for arbitrary a ^ fi. These relations are shown also 
in Appendix In the case of SU(2), {xfx^"') = 
2{siSj). The uniform susceptibility is equal to inverse 
temperature times the Fourier transform of the corre- 
lation function with wave-number zero. Furthermore, 
the specific heat is calculated from the internal energy 

III. HIGH-TEMPERATURE EXPANSION 

In this section, we review high-temperature expansion. 
Although the Heisenberg model is used for explanation, 
this review is more general. That is, one can easily re- 
place the Heisenberg model with a different model. 

A. Moments and cumulants for classical variables 

Fundamental properties of moments and cumulants are 
frequently used in this paper. In order to conveniently 
refer to those properties, let us review moments and cu- 
mulants of classical variables Xi {i = 1, 2, • ■ • ). We put 
most of the fundamental details of the review in Ap- 
pendix 1^, and let us here just show one important for- 
mula used frequently in this paper. Here, a moment is 
denoted by {xjXk ■ ■ ■xi)x, and a cumulant is denoted by 



3 



Moments can be expanded using cumulants, and vice 
versa. However, in those relations, the number of terms 
in the expansion drastically increases as the order of mo- 
ments or cumulants becomes higher. To avoid it, another 
relation, 



{Xi--- Xi) 

-E 



Xi ' ' • Xj\x {-^k ' ' ' •^m) X ; (f ^) 



can be used.0 Here, ^ is the number of a;- variables in the 
bracket in the l.h.s. Then, the summation denoted by 
■p(^, 2) is taken over every partition of ^ elements into 
two groups on the condition that one of the variables, 
for example Xi, must always be included in [•••]^. In 
addition, each bracket includes at least one x-variable. 
In this paper, we call this formula "the mixed expansion 
pivoting on Xi" . For example, the mixed expansion of 
[a;ia;2a;3]a; pivoting on xi is written as 

[xiX2X3]x = (.Tia;2a;3)j; - [a;ia;2]:c(a;3)^ - [a;ia;3]3;(a;2)a; 

-[xi]x{x2X3)x- 



B. Moments and cumulants for quantum variables 

As a reference system, we take a non-interacting sys- 
tem such that Ji_j — for every pair. Its partition 
function is written as Tr 1, and the average in the 

non-interacting system is denoted by 



(O)o 



Tr(6) 
Trl ' 



(16) 



where O is an arbitrary operator. This average (•••)o 
plays a role of (• ■ • )x in the previous subsection. 

There are a number of ways to define moments 
and cumulants of .nuantum variables because of their 
noncommutativity.B One of the ways is defining a mo- 
ment by a symmetrized product. 



Ja{2) 



■ ^V(fc) j<,(fc))o, (17) 



where a represents a permutation of the indeces. Note 
that the ordering of exchange operators Pi j is unim- 
portant in (• • • )s. Hereafter, we use (• • • )o rather than 
(• • • )s if they are obviously equivalent to each other, e.g., 
{Pi.j)s — {Pi,j)o] further examples are commented in Ap- 
pendix 

The partition function, Eq. (10), is rewritten as 



Tr p^/^W 

^ = Trl^;^=Zo(e-^«)c 



(18) 



This averaged Boltzmann factor works as the generating 
function of the symmetrized moments mentioned above. 



namely, with Xi j 

{P,,r ■ ■ Pk,i)s 



d 



9A, 



k,l 



A=0 



The corresponding cumulants are denoted by [• 
given by 



[Pij ■ ■ ■ Pk,l]s 



d 



d 



dX 



t,3 



k.l 



■ln(e-^«)o 



(19) 
, and 

(20) 



A=0 



These moments and cumulants have the properties ex- 
plained in the previous subsection and Appendix 
The free energy is rewritten as 



F ■ 



1 g Ur^rt-u (21) 



^\P^un■■■P^^.3Js■ (22) 



Many terms in the summation of Eq. (^2|) are equivalent 
because the ordering of operators in [• • ■ ]s is unimportant. 
That is, in the case of the cumulant [P-^ j^Pj |^ • ■ • ]s, then, 
m\ I {ki2\k2'i \ ■ ■■) terms are equivalent. 
From Eqs. (|2l| ) and (|2^), we can derive 



OF 

(-/3)"-i 



^ {m-lV. "' , ^ 



' ' ■ Jim-l,jm-l \Pi,jPil,jl ' ' ' Pi,-n~l,ji-n-l\ 



(23) 



m— 



Even if J^j = in the original Hamiltonian {e.g., for 
non-nearest neighbors in the case of nearest-neighbor in- 
teraction), one can start from the general formulation 
written above. After deriving all the formulae, one can 
put Jij- = J or according to the original Hamiltonian. 

Let us consider cumulants for Pij instead of Pi j . See 
the definition of cumulants, Eq. (20). If one replaces P^j- 
by Pi J , the generating function of cumulants changes only 
by — l/?^X](i j) -^ij i ^'^'i thus only the first-order cumu- 
lants are affected by this replacement. In other words. 



differs from 



\P^As - (^j>o - 1/n 



[Pi,i]s = (p»,i)o = 0. 



(24) 



(25) 



On the other hand, the higher order cumulants are equal. 
Namely, 

[Pii-JiPi2,j2 ■ ' ' Pie.Je]s = [Pii.jiPi2,j2 ' ' ' Pie,je]s (26) 
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when [• ■ ■ ]s has two or more pij operators. For practical 
calculation, we mainly use Pij for convenience. In some 
cases, however, the property {pi,j)o — helps simplifica- 
tion. Hereafter, we use Pi j and pi_j interchangeably with 
keeping in mind Eqs. (p^)-(p6|). 



C. Properties of the cumulants 

Let us introduce a diagrammatic representation of the 
moments and the cumulants. Although we mainly use 
those diagrams for later calculation, let us define them 
here in advance because they are convenient also for sim- 
ple explanation in this review. 

Figure |l| shows examples of the diagrams for Pi.j)s 
and [Il-Pijls- We call them bond-diagrams. Each dot 
represents a site. A segment connecting sites i and j 
represents Pi.j, and we call it a bond. Note that a prod- 
uct of two equivalent operators appearing in (■ • • )s [e.g., 
PsaPsa the r.h.s. of Fig. |l](a)] cannot be reduced to 
the identity because these are not always next to each 
other in the symmetrization. Figure |l|(a) is a moment 
bond-diagram, and (b) is a cumulant bond-diagram. 



(b) 



[-Pl,4-Pl,2-P2,3f3,4]s 



FIG. 1: Examples of bond-diagrams. 



A cumulant is equal to zero if the variables in [■ ■ ■]s 
can be partitioned into two groups which are indepen- 
dent of each other on averaging (•••)o- For example, 
[Pi_2^3,4]s = because the trace of sites 1,2 is taken in- 
dependently of that of sites 3,4 in the non-interacting 
system. In other words, a cumulant is equal to zero if 
the bond-diagram is not connected as shown in Fig. ^(a). 
Here, a solid rectangle in Fig. |^ represents arbitrary 
bonds. Therefore, as for cumulants, we consider only 
connected diagrams hereafter. 



(a) 



= 



(b) [r-,]--=[Pi,2]s = {Pi,2)o = - 



n 



(c) 



= 




FIG. 2: Properties of the cumulants. A solid rectangle rep- 
resents arbitrary bonds. 

Figure ^(b) is a first-order cumulant. In a higher-order 



cumulant bond-diagram, more than one bonds must meet 
at each dot. If not, it is equal to zero as shown in Fig. ||(c) 
and (d). Furthermore, a cumulant bond-diagram is equal 
to zero if two parts in it are connected by only one bond 
as shown in Fig. |^(e). We prove these properties in Ap- 
pendix 



D. An example of the finite cluster method 

The finite cluster method (FCM) is [a-atandard method 
for the high-temperature expansion.Eia It is also used 
for perturbajtive expansions starting from exactly solv- 
able modelsJj In order to explain it in detail, we give 
an explicit example using the one-dimensional nearest- 
neighbor interacting system with a periodic boundary 
condition, namely, Jij — J(5j.i+i, (1 < « < N), where 
the site iV -I- 1 is equivalent to site 1. In this case, the 
procedure of the FCM is reduced to a very simple form. 

In the FCM, series coefficients in the thermodynamic 
limit are exactly obtained by summing contributions 
from finite-size clusters. The clusters are subsets of 
the infinite lattice, and we need only connected clusters 
because cumulant bond-diagrams must be "connected" . 
Namely, the one dimensional system requires only open 
chains. The cluster of size i is defined as 



H 



(p) 



/ 1\ 



(27) 



where we have chosen to use pij rather than Pij in order 
to make use of [pij]s = 0. The expectation value in this 
cluster is given by 



{Pi,j)i ■■= 



') 



E 



Pt,j {n 



(28) 



We assume i < j hereafter. Note that {pij)i = when 
i < 1 OY j > I because of a property shown in Fig. ||(c) 
and (d). Since a finite cluster is used instead of the origi- 
nal system, {pij)e contains only a subset of terms appear- 
ing in {pi,j). The important point is that {Pi,j)£ includes 
contributions from clusters smaller than £, and thus sim- 
ple summation of (jpi.j)i over I yields multiple counting 
of those terms. To avoid this, the FCM requires subtrac- 
tion of contributions from the smaller clusters. We define 
a net contribution from the £-site cluster, 

£1=0^2=0 

where the summation excludes (-i = (.2 = 0. Here, 
(■I (^2) represents the number of removed sites from the 
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side of 1 {€). When i = \ and z + a; = this expression 
is reduced to (pi.f)^ ~ {pi.e)e- Let {pi.i+x)<e denote the 
total net contribution to (pi^i^^) from clusters smaller 
than or equal to £, 



{Pi,i+x)'< 



Y{Pj:3+x)'l- 

l=x+l j=l 



(30) 



Here, {pi^i+x)'<e does not depend on i because of the 
translational symmetry. Note that (pi.i+x) — {Pi,i+x)'<oo- 
When £ is finite, the series for {pi^i+x)'<g is correct up to 
a finite order. To think about the "finite order", let us 
take a cumulant bond-diagram appearing in the m-th or- 
der. The diagram has a bond Pis+x and m bonds from 
the Hamiltonian. In order to give a nonzero contribution, 
there must be at least one bond per nearest-neighbor pair 
between i and j. Furthermore, nearest-neighbor pairs 

not between i and i + x less than i or greater than 

i -\r X must have at least two bonds per pair. There- 
fore, there is a contribution from l)-site cluster when 
m > X + 2(£ — x). That is, {pi^i+x)'<_i is correct up to 
0[(/3J)2^-"-i]. 

The above is the usual method of the FCM. However, 
in one dimension, these equations can be reduced to sim- 
pler forms as shown in Appendix 



IV. A NEW METHOD FOR THE 
CORRELATION FUNCTION 

Let us think about calculating the Fourier transform 
of the correlation function, 



(31) 



Series of S^{q) up to 0[{PJ)^] requires series of each 
{pj+x,j) np to 0[{PJ)^]. Here, the important point is 
that the series coefficients of {pj+x,j) up to 0[(/3 J)'"^'"^] 
are equal to zero as explained in Fig. ||(c),(d) and in 
Sec. [HD. Hence, we need (pj+xj) only for 1 < a; < M. 
In addition, that property provides us a weak point of the 
FCM: When M is fixed, large x requires a larger cluster 
than small x does. Thus, the most time-consuming part 
is calculation in the largest cluster, x = A/, namely, cal- 
culation of Tr(i+^){pi,i+M(H[^j\^)*^}. However, what 
we need for x = M is only the M-th order coeflScient be- 
cause we already know that the lower orders are equal to 
zero. Therefore, if only the FCM is used, the most time- 
consuming calculation yields very little information. 

Our goal in this paper to formulate a new method is 
as follows. If the contributions from large clusters can 
be calculated by anothet. method, the FCM can be used 
only for smaller clusters .Q Accordingly, one can calculate 
up to high orders. In Ref. ^, we have calculated S^{q) 
up to 0[(/3J)^^]. However, we have used the FCM only 
for i < 13. The corrections to the results have been 
calculated by the method explained in this section. 



A. Traces using combinatorics 

We calculate traces of products of exchange operators 
by decomposing every permutation into a product of in- 
dependent cyclic permutationsDEl as explained in the fol- 
lowing, (an explicit example is given in Appendix |^. Let 
us consider a trace of P := Pii,jiPi2,j2 ' ' ' Pimdm^ with 
Tr*-^-* denoting the trace in the ^-site system, 



n n 

/ 1 2 f I I 1 2 I 

y. ' ' ' / . \Q!iQ^2 ■ ■ ■ ai \-r\ ctxoii ■ ■ ■ at) 



ai — 1 a^ — 1 



■ ■ ■ \pLi(X2 ■■ ■ ai\ apiap2 ■ ■ ■ ape) 

ai — 1 a^ — l 
n n 

Ql — 1 Q£ — 1 

where api refers to the color at position i after the per- 
mutation P. The summation of ai makes a contribution 
only when ai — api for every i. Consider using this re- 
lation successively starting from i. That is, ai is equal 
to api, and then api is equal to apsj, . . . , one can re- 
peat this procedure until coming back to ai at a certain 
power of P, namely, a; — api — ap2j — ■ ■ ■ — on. In 
other words, all the variables whose subscript belong to 
one cyclic permutation in P have to be equal. Since any 
permutation can be decomposed into a product of in- 
dependent cyclic permutations, the number of indepen- 
dent variables of the summation is the number Y{P) of 
cyclic permutations of P. Therefore the trace is given by, 
p = n^(P\ and accordingly. 



{P)i = 



TrWp 
Tr^l 



(33) 



In fact, this method allows us to calculate Tr^^\Pij P 
at the same time as Tr*^^'P for an arbitrary pair.l 
Let us remember the permutation P is decomposed into 
cycles in calculating Tr^^-'P. When the cycle which i be- 
longs to is different from what j belongs to, the next oper- 
ation Pij unites the two cycles into one. Then, the num- 
ber of cycles decreases by one, namely, Tr'^^ (P^^P) = 
nX(P)-'^, On the contrary, if i and j belong to one cy- 
cle of P, then Pij breaks the cycle into two. Then, the 
number of cycles increases by one, namely, Tr^^-'(Pi jP) = 



B. Explicit calculation of cumulants 

To represent {YlPi,j)o: we introduce an unsym- 
metrized version of a bond-diagram. Figure ^ shows its 
examples. Since this diagram is equivalent to an "amida 
lottery", used in Japan to peacefully decide how to dis- 
tribute a fixed number of objects among an equal number 
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of people, here we call it an amida-diagram. Each vertical 
line represents a site, and a horizontal line between site i 
and j represents Pij. The order of horizontal lines from 
bottom to top should coincide with the order of Pij in 
(• • • )o from right to left. We bracket an amida-diagram 
by (• • • ) in order to distinguish it from a cumulant version 
of the amida-diagram introduced later in Sec. 



(a) 



(b) 



1 Z .5 ^ 



{-P3,4-Pl,2f3,4f2,3fl,2-P2,3)o 



FIG. 3: Examples of amida-diagrams. 



(a) 



(b) 



n 

1 ( 



■1) 



n 



FIG. 5: The mixed expansion pivoting on Pi_f for (a) the 
lowest-order and (b) the second-lowest-order nonzero cumu- 
lants in the (pi,f ) series. The dotted line represents a sequence 
of bonds such that there is only one bond between sites. 



Since expectation values are not changed by transposi- 
tions of neighboring commutative operators, the ordering 
of Pi, 2 and P3,4 is irrelevant to the results. On the other 
hand, two exchange operators sharing one of the site- 
indeces, such as Pi, 2 and P2,3, are not commutative, and 
in general the ordering of those bonds affects the result. 
For example, the value of Fig. ||(a) is different from that 
of (b). However, in some cases, even the order of such 
"site-sharing" bonds does not affect the result. Let us 
think about the value of an amida-diagram in Fig. |^(a) 
by counting cycles. The "initial" state is the identity 
permutation, which has H cycles. Then, every time Pij 
is applied, it unites two cycles into one. After all the 
Pi.j operators are applied, there is only one cycle, which 
includes all the £ sites. Hence the expectation value is 
n}~^ , and it does not depend on the ordering of those 
Pi,j operators. Therefore, after the symmetrization, the 
moment bond-diagram Fig. ^(b) also gives n^~^. 

(w (—■■.■■—>-»'"' 

FIG. 4: Actual calculated values of diagrams. The dotted 
line represents a sequence of bonds such that there is only one 
bond between sites. 

Let us think about adding Pi^ to Fig. ^(a) on the top. 
Since sites 1 and £ belong to the same cycle, the oper- 
ation of Pi,^ cuts the cycle. Namely, the permutation 
has two cycles and the diagram gives n^~^. Again, af- 
ter the symmetrization, the corresponding moment bond- 
diagram also gives n^~^. 

In Fig. ^(a), we calculate a cumulant using the values 
derived above. This is the lowest-order nonzero cumu- 
lant in the {pi/) series. We use the mixed expansion, 
Eq. ([isl), pivoting on Pi/. The expression of the expan- 
sion is very simple because [Pi,^]s is the one and only 
nonzero cumulant which includes Pi/. 

Next, let us consider adding one more bond. Pi, 2- The 
mixed expansion pivoting on Pi,^ is shown in Fig. |^(b). 
The coefficient 2 of the last term is due to the two ways 



of choosing which one of the two Pi, 2 is put in [■ • ■]s. 
The first term of the r.h.s. is calculated as shown in 
Fig. ^ The symmetrization is equivalent to averaging 
over all the possible configurations of bonds. In fact, 
however, the value of each amida-diagram is determined 
only by relative configurations of two Pi, 2 and one P2.3. 
There are only two patterns, namely. Fig. |^(b) and (c). 
Their r.h.s. are obtained by simplification Pi,2Pi,2 — 1 
and Pi,2P2,3Pi,2 — Pi, 3- Here, an "f]-shape" bond in 
Fig. ||(c) represents that the bond is not connected to 
the site at "O". We already know how to calculate the 
r.h.s. of Fig. ^(b) and (c). We also know that configu- 
ration of Pi,i+i [i > 3) does not affect the value of the 
amida-diagram. That is, the value of the amida-diagram 
is determined only by the position of P2 3 . 

In order to obtain the l.h.s. of Fig. ^(a), we have to 
count how often each of the configurations Fig. ||(b) and 
(c) appears in the symmetrization. For example, the pref- 
actor 1/3 in the last term of Fig. ^(a) is obtained by 

{the number of configurations of Fig. ^(c)} ^^^^ 
{the total number of possible configurations} 

It can be considered like this: The vertical line of site 
2 is partitioned into three regions by the two Pi. 2. The 
probability of finding P2,3 in the intermediate region is 
1/3. More explicitly, it can be given also by integrals, 

/odria/' dT{2r-dr23 1 
/o dri2 /^^^^ dr{2 /o dT23 3' 

where T12, t{2, T23 represent positions of lower Pi, 2, upper 
Pi,2, P2,3, respectively. Here, t = is the bottom, t — 1 
is the top. On the other hand, the "probability" of the 
other configuration Fig. ^(b) is 2/3. Consequently, the 
symmetrized value is calculated as shown in Fig. ra(a). 
This formula is true only when £ > 3. li £ — 2, we cannot 
take an average over position of P2,3, and Fig. ^(c) does 
not appear. 

Examples of higher order results are shown in Fig. |^. 
We can find simple properties in these cumulants. 

(i) A cumulant is reduced by a factor l/n as £ increases 
by one. 
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FIG. 6: Calculation of a symmetrized moment using amida- 
diagrams. 



Pk^i in the l.h.s. unites cycles or breaks a cycle, Pj,i in 
the r.h.s. does the same thing. Therefore, the number of 
cycles of the l.h.s. is equal to that of the r.h.s. As for 
the denominator of the trace, the l.h.s. has one more site 
than the r.h.s. Hence, in total, the r.h.s. needs a factor 
l/n. The logic above is true even if site i is the same as 
site I. 




(b) 
(c) 
(d) 



1 2 f 



- 1 /n 2 



1 2 e 



1 2 t 



12 t 



\ 3 n 



n'-l 6 
—\-^+ — 



i« \ 30 













n J 



r?- \ ( 1 



,f \ 3 

FIG. 7: Actual calculated values of cumulants 



{I > 3) 
{£ > 3) 
{£ > 4) 
> 5) 
{£ > 2) 



FIG. 8: A part of an amida-diagram. A bond can be deleted 
to yield a factor l/n. 

The symmetrization of all the operators except Pj ^. 
is equivalent to the symmetrization of all the operators 
as explained in Appendix Hence, in all the terms in 
the symmetrization, one can put Pj^k at the bottom of 
amida-diagrams and apply the contraction rule above. 
Therefore, we can obtain a contraction property also for 
a moment bond-diagram as shown in Fig. O. 



■ arbitrary 



(ii) A double bond in Fig. ^(a) yields a factor of — ^), 
two distant double bonds in Fig. 0(d) yield square 
of this factor. In fact, we can make a more general 
statement: When multiple bonds are distant, we 
can simply make a product of contributions from 
each of multiple bonds. 

We more rigorously derive property (i) in the next sub- 
section and property (ii) in the subsection after the next. 
After that, we formulate the new method using those 
properties. 



FIG. 9: Contraction of a moment bond-diagram. 



else 



(Theorem I) 



C. The "contractible" property 

First, we consider a contraction of an amida-diagram. 
In an amida-diagram, let us focus on a part that includes 
two sites as shown in Fig. ||. Because of a property of a 
trace, the ordering of operators can be cyclically rotated. 
In other words, amida-diagrams have a periodic bound- 
ary condition. Hence we can put a bond Pj^k at the 
bottom of the amida-diagram without loss of generality. 
Let us compare the l.h.s. and the r.h.s. at each position 
of bonds. In the l.h.s., the bottommost bond Pj^k unites 
two cycles at sites j and k of the initial state, and makes 
a cycle of the two sites. At this stage, nothing happens 
to the r.h.s., and there is a cycle at j. Then, the rest of 
the operations of the l.h.s. and the r.h.s. are the same. 
That is, Pij unites this cycle and another cycle; when 



FIG. 10: Theorem I and its proof. A solid rectangle represents 
arbitrary bonds. In the expansion, the bonds are partitioned 
into two groups, which are represented by split rectangles. 

Finally, we prove that a cumulant bond-diagram has 
the same contraction property as the corresponding mo- 
ment bond-diagram does. We use a mathematical induc- 
tion for the proof. 

(i) The simplest nonzero cumulant, the l.h.s. of 
Fig. |^(a), has this contraction property. 

(ii) Let us think about a certain cumulant bond- 
diagram. We assume that "the lower-order cu- 
mulant bond-diagrams have this contraction prop- 
erty." Then, we use the mixed expansion, Eq. ([T^), 
pivoting on one of the bonds as shown Fig. 11^. We 
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use this assumption for the third term. The prop- 
erty of moments can be apphed to the first and 
second terms. 

Because of (i) and (ii) , the contraction property is true for 
any cumulant bond diagram that satisfies the condition. 
We call this relation the Theorem I. 



D. The "detachable" property 

Some of the cumulants can be decomposed into a 
"background" and local quantities. We give an exam- 
ple in Fig. |ll|. 



( — )q := 



(■ 



2l] X (^)q X 



1. Proof of a lemma 

For the proof, we introduce some notation. It is 
schematically shown in Fig. O. A diagram is denoted 
by G. Here, solid rectangle AJB) represents "rectangle 
I (II) and a sequence of single bonds" in Fig. |l^ unit- 
edly. We define Ga (Gb) as a diagram such that each 
multiple bond in B (A) of G is replaced with a single 
bond. We define the "background" Ag as a cumulant 
bond-diagram such that each multiple bond in A and B 
of G is replaced with a single bond. Hence, there is a 
relation, Ag Ag^ = ■ 



Ga 



An 



Gb 



FIG. 11: A cumulant can be decomposed into a "background" 
and local quantities. This property is generally described by 
Theorem II below. 

As already appeared several times, an intersite that 
has only one bond plays a special role in this paper. We 
call it a single bond here. On the other hand, an intersite 
that has more than one bond is called a multiple bond. 

In this subsection, we prove Theorem II shown in 
Fig. ^ Here, a solid rectangle represents arbitrary 
bonds. A dotted line stands for a sequence of single 
bonds. The bond-diagram in the r.h.s. of (a) has the same 
length as that of the l.h.s.; the bond-diagram in a denomi- 
nator of (b) has the same length as that of the numerator. 
This decomposition is possible even if the bond-diagram 
has more than two blocks of multiple bonds, on condition 
that the number of single bonds should be greater than 
that of blocks. 



(a) 



(b) 







c 




m»mi 




1 


'-I 











FIG. 12: Theorem II. The details are explained in the text. 

In other words, the Theorem II is described as fol- 
lows: a bond-diagram can be rewritten by the product 
of a background and local quantities. The background is 
obtained by replacing each multiple bond with a single 
bond. One of the local quantities is obtained by replac- 
ing each multiple bond with a single bond except those 
in that block, and dividing it by the background. 



FIG. 13: Definitions. A dotted line represents a sequence of 
single bonds. The details are explained in the text. 



Let us consider the cumulant of G - 



As 

is also shown in Fig. 14, the mixed expansion pivoting on 
is written as. 



[G], = Ci(G) + C2(G), 



Gi(G) := {G). - [P^,,]s{l[Pk,i 

k,l 



G2(G) 



else k' ,1' 



Pk' 



k'V 



k".l" 



(36) 



(37) 



(38) 



That is, the second term of Ci is decomposition into Pij 
and the rest. The other decompositions are included in 
G2. In G2, each cumulant includes more than one ex- 
change operators, and thus, to give a nonzero cumulant, 
each site index has to appear at least twice as explained 
in Fig. |^(c) and (d). In other words, each cumulant bond- 
diagram in C2 has to have a loop. 



^Ci(G)+C2(G) 



Ci(G) 
C2{G) 



FIG. 14: The mixed expansion pivoting on Pij . The bonds 
are partitioned into two groups, which are represented by split 
rectangles. 
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In order to prove the Theorem II, we first prove a 
lemma, 



Ci(G) Ci(G^)Ci(Gb) 



Af 



(39) 



in the following. 

First, we give an explicit examp le of Ci using the cu- 
mulant in Fig. |^(b). In Sec. IV B , we have symmetrized 
each of the terms in the mixed expansion individually. 
Here, however, we postpone the symmetrization proce- 
dure. That is, as shown in Fig. 15, we first sum amida- 



That is, as shown in Fig. 
diagrams of the terms in Ci with a certain configuration 
of bonds. After that, we symmetrize by averaging over 
all possible configuration of bonds. We carry out the 

symmetrization of all the operators except Pi^i this 

is equivalent to the symmetrization of all the operators 
as shown in Appendix |c[ Let us consider Fig. ^(a). In 
the amida-diagram in the second term, sites 1 and £ be- 
long to different cycles. Hence, when one more bond Pi ^ 
is added (to make the amida-diagram of the first term), 
the two cycles are united into one. Thus, the ratio of 
these amida-diagrams is 1/n, which is equal to [-Pl,^]s in 
the second term. Accordingly, these terms cancel out 
each other. Therefore, there is a nonzero contribution to 
Ci only when two sites 1 and £ belong to one cycle in 
the "amida-diagram without Pi/" in the second term of 
Ci. Hereafter, such a configuration is called a Ci-finite 
configuration. Figure |l5|(b) is a Ci-finite configuration. 
Its contribution to Gi (Gpig E|(b) ) is determined by "devi- 
ation from the background . There is an "extra" bond 
compared to the background, and here it increases the 
number of cycles by one to make a factor n. On the 
other hand, the "probability" for this configuration to 
occur is 1/3 as shown in Eqs. (^J) and (35). Therefore 
we obtain Gi(GFig.|b)) = ^nAa^^^i,)- 



(a) 



(b) 




FIG. 15: 

The example above is a simple case. In general, "the 
probability of the configuration" times "the deviation 
from the background" is summed over all the Gi-finite 
configurations, namely, 

Gi(G) = Ag ^ (probability) x (deviation). (40) 

Ci -finite 
config. 



Let us consider a diagram G shown in Fig. |l^. Then, 
in order to put sites 1 and £ in one cycle in the "amida- 
diagram without Pi/" of G, each of Ga and Gb must 
also have a Gi-flnite configuration. Then, the deviation 
from the background as a whole is given by (deviation in 
A) X (deviation in B). 

In fact, "probability" is also a product of those of A 
and B as shown below. For example, the probability of 
configuration in Fig. |l^(a) is written as 



2" / dTi2 



<^Tl2 / dT23 

Tl2 •^'ri2 

dT45 / dr^s / dT34. (41) 

Jr45 J Tis 



The first three integrals are independent of the last three 
integrals, and thus the probability is written by the prod- 
uct of two parts. Next, let us think about Fig. [lq(b). 
If we naively use r, the expression of the integrals is 
complicated. However, if we utilize the periodic bound- 
ary condition of r and shift the integration variables as 
p := T - T23 for r = r34, t^^, we obtain. 



dTi2 / dr^j / dr23 
J T12 J T12 

d/334 / dp34 / d/945- 

J P34 J P34 



(42) 



This expression is equivalent to Eq. ( p4] ) . Hence the prob- 
ability is written by the product of two parts again. The 
important point in deriving Eq. ( p2| ) is existence of P4,5. 
Suppose the diagram does not have P4,5, and has Pi_4 
instead of Pi, 5. The uppermost bond Pi, 4 is always at 
T = 1. If the integration variables are shifted, 1 — r23 
appears in lower or upper bounds of the integrals. That 
is, the integrals are not written by a product of two parts 
any more in that case. 



1 2 3 4 5 



(a) 



(b) 



1 2 3 4 5 



FIG. 16: 



The shift of the integration variables mentioned above 
can generally be applied to any blocks A and B. Hence, 
the probability in Eq. (^ is also a product of those of 
A and B. The summation over Gi-finite configuration is 
double summation over Gi-finite configuration of A and 
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over Ci-finite configuration of B. Tlierefore, 



Ci(G) 
Ag 



(probability^) (deviation^) 

nite 

s.A 

X ^ (probability i?) (deviations) 

Ci-tinito 
contig.B 

C^{Ga)Ci{Gb) 



Ci-finite 
contig.A 



^Ga 



(43) 



Acg , we obtain the lemma, 



By usmg Ag = Ag^ 

Note that the shift of the integration variables men- 
tioned above can be applied even if a diagram has more 
than two blocks of multiple bonds. A block and a sin- 
gle bond is regarded as one group whose members have 
a common shift of the integration variables. This is the 
reason for the restriction of the Theorem II; the number 
of single bonds should be greater than that of blocks. 



2. Proof of the Theorem II 

In order to focus on the "multiple bond" region, we 
introduce another notation for moments and cumulants. 
Let us think about Ga- A set of bonds. A, is partitioned 
into two subsets in calculating C2{Ga)- Let a denote a 
subset of A. Then, Q{a) denotes the cumulant in which 
A is replaced by a. In choosing a, at least one bond 
per intersite must be chosen to make Q{a) nonzero. A 
moment for bonds A — a is denoted hy Ai{A — a). Then, 
Eq. (p6|) for Ga can be rewritten as 



Q{A) = CiiGA) 



Q{a)M{A-a). 



(44) 



The l.h.s. is equal to the excluded term a = Ain the sum- 
mation in the r.h.s. because A4{9) = (l)o = 1- Therefore, 
we obtain 



Ci{Ga) - ^Q(a)A4(A-a) 



(45) 



The corresponding equation for graph G is 
Ci{G) ^ 



A 



+ B-b). (46) 



(47) 



Since A is not connected to B, 

M{A -a + B-b)^ M{A - a)M{B - h). 
On the other hand, according to the lemma and Eq. (45) 
Ci(G) Gi{Ga)Gi{Gb) 



EE 



A 

Q{a) Q{b) 



M{A~a)M{B -b). (48) 



Here, we use a mathematical induction again for the 
proof. 

(i) The Theorem II is true when there are two extra 
bonds on the background as explicitly calculated in 
Fig. 0(d). 

(ii) We assume that the Theorem II is true when the 
number of extra bonds on the background is less 
than that of G. According to the assumption, 



Q{a + b) Q{a)Q{b) 



{a + b^A + B) 



(49) 



Combine the r.h.s. of both Eqs. (46) and (f48|), use 
Eqs. (^) and (|9|). Then, all the terms except for 
(a, b) = {A, B) cancel out and we obtain 



Q{A + B) ^ Q{A) Q{B) 
A ~ A A ■ 



(50) 



According to (i) and (ii), the Theorem II is true for any 
G that satisfies the conditions. 

The proof mentioned above can be applied even if block 
AoT B includes single bonds. Then, the Theorem II for a 
diagram with two blocks is repeatedly used to prove the 
Theorem II for a diagram with more than two blocks. 
Therefore, also a diagram with more than two blocks has 
the "detachable" property as mentioned first. 



FIG. 17: Theorem F 



3. Theorem F 



Now that the Theorem II is proved, the contraction 
can be written in a more general form. Although the 
Theorem I is for a sequence of three single bonds, se- 
quences of two single bonds also can be contracted as 
shown in Fig. ^ as Theorem F, on condition that each 
block must be separated by at least one single bond and 
two of the blocks must be separated by a sequence of at 
least two single bonds. The Theorem F is easily proved 
using Fig. |2|(a); contract Ag in the r.h.s. and reversely 
use the same relation to the contracted diagram. 



E. The quasi-cumulant method 

Apart from the restriction, the Theorem II is equiva- 
lent to a property of moments a disconnected mo- 
ment is equal to a product of connected moments. In 
order to effectively use the Theorem II, we define quasi- 
moments and corresponding quasi- cumulants below. Our 
definition of a quasi-moment is a factor times a "local 



11 



quantity" in the Theorem II. Although another defini- 
tion can be the local quantity itself, it is less convenient 
for our purpose as discussed in Appendix |^. 



(-Pl,2>q 



1 [^] 



2 [^] 

1 [^^] 



(Pi 



(o)^2i 



3-2 [rrr^] 



2 [ 

FIG. 18: Examples of quasi-moments. 



Only up to a finite order, the Theorem II can be applied 
and Qqm works as a generating function. The larger 
is, the higher order the Theorem II is valid up to. Here, 
we assume that |i — j| is large enough, and forget about 
this restriction for the time being. It will be commented 
on later. Then, the quasi-moments are defined as 



k,l 



■ Py,^)^ 



dX 



k,l 



dXy 



(52) 



A=0 



Let us think about the one-dimensional system again. 
We define the generating function of quasi-moments as 



"qm 



1 



[PiJ Pi,i+lPi+l 



1+2 ■ 



■P,-l. 



(51) 



The zeroth-order quasi-moment satisfies (l)q = 1. Some 
examples of quasi-moments are shown in Fig.|l^. More 
explicitly, the quasi-moments are obtained by 



-^1,2 ^2,3 



pmi + 1 pm2 + l 
^1.2 ^2.3 



P, 



,mfc-l+l 



k-l,k 



Pi 



k,k+lPl,k+l 



(53) 



q (rrii -M)(m2 + 1) • • • (rrifc-i + 1) [Pi,2-P2,3 • ■ • -Pfc-i,fc-Pfc,fc+iA,fc+i]s ' 
where we have assumed that all the bonds are on the "background" . We need a special care when Pk i is outside the 



background, namely, /c < i or j < Z in the case of 5qm • Such a operator is explicitly denoted by P^°) . Then, we can 



more generally write as. 



(p(o))™,...(^{o) ) 



'Pi 



■P, 



k-l,k 



\-'^1.2 ) ' ' ' 



(o) 



p™fc-i + lp p 
Pk^l.k t^k,k+in,k+l 



(TOi -t- 1) • ■ • [nik-i + 1) [Pis+i ■ ■ ■ Pk-i.kPk,k+iPi,k+i] 
I 



(54) 



Next, using quasi-moment defined above, the generat- 
ing function of quasi-cumulants is defined as. 



:= InC**'-''). 

^qc ^qm 

Then, the quasi-cumulants are defined by, 



[Pk,l ■ ■ ■ Py,z]q ■ — 



d 



dX 



k,l 



(55) 



(56) 



A=0 



According to the Theorem II, when none of fc or Z is equal 
to y or z. 



exp 



exp 



Xk,lPk,l + ^ Xy^zPy.z 



k,l 



Xk,iPk,i 



kl 



exp 



(57) 



and thus, as discussed in Appendix 



k,l 



= 0. 



(58) 



Therefore, only "connected" quasi-cumulants remain 
nonzero. If only the quasi-moments are used, we have 
to classify connected and disconnected quasi-moments. 
However, if the quasi-cumulants are used, we need only 
connected quasi-cumulants and can avoid the complicated 
classification. 

Suppose we set A^^; — X for every (fc, I) pair. Then, 
g(^^i+x) ^ 0) is a function of x and A. The 
quasi-cumulants have a translational symmetry in the 
background. For example, [P™+i]q = [f^+i^^+2]q = 
■ ' ' — [Pi+x-i i+x]q' ^^'^ thus, these "length-1" quasi- 
cumulants can be summed up to a; [-Pj™+i]q in Sqc*"'"^'*. 
If we in the same way simplify Gqc^~^^\ length-2 quasi- 
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cumulants are multiplied by x — 1, length-S quasi- 
cumulants are multiplied by x — 2, and so on. The quasi- 
cumulants including pjfi are located at both ends, and 
thus, they are multiplied by two. Namely, 



(59) 



l<ml<m 



+(■^-1) E 



1-^7,2+1-' 



l<mi ,m2 
mi +m2 <m 



+ 

+ {x-m+ l)A"[Pi,,+i • ■ • P,+m-l,z+m]q, (60) 

2<mi <m 

^ , \mi+m2 

-2 E ^[(i^^l.)-^"ir]. 



2<mi,l<m2 
mi+m2 <m 



\mi+m2 



2<mi ,m2 
mi +m2 <m 



mi!m2! 



[(/^r2..-l)""m.J""]q 

(61) 



Hence, this expression for arbitrary x can be derived if we 
obtain all the quasi-cumulants up to m-th order, which 
are calculable in small systems. Therefore, we can obtain 
{Pi,i+x) up to (a; + TO)-th order for arbitrary x by tracing 
back Eqs. (||)-(|ll), (H) and We call this method 

the quasi-cumulant method (QCM). 

The QCM has been done by a brute-force program 
of Mathematica. The quasi-cumulants can be calcu- 
lated by Eq. (|5^), that is, differentiating the series of 
their generating function with respect to corresponding 
Ji^i+i, where Ji^i+i are coupling constants of the Hamil- 
tonian of a rather general form, Ji,i+iPi,i+i- We 
take the system as small as possible on condition that 
it be large enough to correctly calculate those quasi- 
cumulants. That is, quasi-cumulants of length I can be 
obtained in the system of length 1. First, we calculate 
the series of {pij) to obtain ^/qm'' up to m-th order. At 
each order of calculating Ti.'^ , we exclude terms irrelevant 
to the derivative in order to save computational time and 
memory. Let us rewrite 5qm •* as 



k=0 



k\ 



(62) 



then, Eq. (|T5|), the mixed expansion pivoting on one of 
Ti. can be used, 



k-l 



(fc-l)! 



-[{i-mk-iy. 



[H']q(H'=-')q, (63) 



E 

k=l 



kl 



0(/3'"+i). (64) 



Now we are ready for the derivatives. The point is that 
we need at most m+2 sites here, while we can reconstruct 
the generating function of larger systems from the quasi- 
cumulants and thus we can obtain {pi^i+x) up to {x + m)- 
th order arbitrary x. 

Finally, we comment on the restriction of the QCM. 
The QCM is based on the Theorem II. Let us consider 
diagrams shown in Fig. We cannot apply the Theo- 
rem II to these diagrams, and thus Eq. ( |57| ) is not valid. 
However, we can apply the Theorem II if the order is 
lower than these. Hence, we can use Eqs. (^9|)-(|6l|), only 
for m < x/2 {x: even) or m < (a; — l)/2 [x: odd). 



(b) (^^^ir>^ 



FIG. 19: The Theorem II is not valid for these diagrams. 

Note that the QCM is useful for contributions from 
large clusters rather than those from small clusters. 
Hence, by combining the QCM with the FCM, we can 
calculate the high-temperature series coefficients up to a 
high order. Figure ^ shows which of the FCM and the 
QCM was used to obtain the series coefficients in Ref. IJ. 
We calculated {pi^i+x) up to (a; -I- 6)-th order for x > 13 
by the QCM. Note that the QCM results are vahd for 
arbitrary x > 13, including the x ^ oo limit as indicated 
in Fig. The reason why we stopped the QCM at this 
order was mainly due to the FCM. If the FCM can treat 
larger systems, the restriction of the Theorem II is less 
important, and the QCM can go further. Our FCM pro- 
gram was for general SU(n); we expect that more terms 
can be obtained if the interest is restricted to SU(2) in 
making the FCM program. 

In Ref. ^, the largest cluster we used for the FCM was 
composed of ^ = 13 sites. As noted in the context of 
Eq. @), {p^,^+x)'<i is correct only up to 0[(/3J)2^-"-i], 
and higher orders include contributions from the larger 
clusters. Hence, for high-order coefficients of x < 12, 
corrections to the FCM results must be calculated. The 
QCM can be used also for this purpose because these cor- 
rections are a subset of the terms obtained by the QCM. 
We put a fictitious coefficient on each quasi-cumulant 
located at both ends in Eq. (|6ll). These fictitious co- 
efficients carry information about the cluster-size of a 
cumulant, and we can distinguish terms needed for the 
corrections to the FCM result. Coefficients obtained by 
this procedure in Ref. | are indicated by "FCM+QCM" 
in Fig. 1^. 

Since these series coefficients are obtained as explicit 
functions of n, the data of the series coefficients are too 



13 



{Pi,i+x) distance x 

1 2 1213 19 



1 

2 

order 








FCM h 








QCM 






numerous to be fully listed in this paper. Hence, we list 
here only the results of the QCM, (pj^j+j,) up to {x + 6)- 
th order, and the other data will be provided on request. 
The coefiicients of the series 



{Pi,i+.) = n-'-- (-1 + n^) n,,m (-/3 J)™ (65) 



are given in the following: 



1 {x > 1), 



(66) 



FIG. 20: We obtained high-temperature series coefficients of 
{pi,i+x) using the FCM and the QCM complement arily. The 
figure shows which motliod was used for oacli scries coofhcicrit. 
Though we have used only all the coefficients up to 0[(/3J)^®], 
the QCM results are valid for arbitrary x > 13. 



(-6 + n'^) X , , , ^ 

o n 



360 



5a:2(-6 + n2) -6a; (-30 + 15 + 2 n^) + (-150 + 7n^ 

[x > 4) 



(68) 



TTxx+s = — 2,bx^ i-& + n^f -U&x^ (l80- 120n^ + 3n^ + 2n'') 

' 45360 L ^ ' ^ ' 

+2 X (-7560 + 17010 - 2961 + 302 n^) 
+24n2 (945 + 2171^-20 n^)l (a; > 6), 



(69) 



n^,x+A = ^^^^ [l75a;^(-6 + n2)''- 1260a;3(-6 + n2)^ (-30 + 15 71^ + 2n*) 

-|-2a;^ (1247400- 1965600 n^-l- 519750 n*- 48210 n^-l- 6817 n^) 
-6 X (-226800 + 1020600 - 297990 - 36795 + 5468 n^) 

+9 (-352800 + 212520 - 14850 + 3287 n^) ] {x>8), 



(70) 



359251200 



385 .x-"^ (-6 



4620 



n^f (-30 



15 + 2n^ 



+44 x^ (-2381400 + 3572100 - 1260630 n"^ 



147915 - 16461 + 2021 n^°) 
-66 x^ (2268000 - 6350400 ri^ + 2891700 - 101940 n*^ - 53013 + 6496 n^°) 
+2 a; ( - 35925120 + 264448800 - 182037240 n'^ 
+ 17223030 n^' - 3830673 + 517472 n^°) 

-36 (-6652800 + 7983360 - 60720 - 165099 + 27776 n*) 
(a; > 10), 



(71) 
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n-6 



TT 



x,x-\-6 



5884534656000 
-3153150 (-6 + (-30+ ISn^ + 2n^) 

+5005 (-6 + n^) ^ (3855600 - 4725000 + 982800 n'^ + 18870 + 18947 n®) 
-6006 ( - 306180000 + 731430000 - 413229600 + 58660200 n*^ 

+2763810 71** - 1239825 71^° + 126908 n^'^) 
+13 (172260950400 - 745940210400 + 577975305600 - 101282227200 

+6937664580 - 2448983250 + 265586033 n^^) 
-6x{- 163459296000 + 1777619844000 - 2058679022400 

+396577981800 + 31224022830 7i^ - 11642353905 + 1392291386 n^^) 
+ 135 ( - 32691859200+ 61509127680 - 15679984320 69539870471*^ 

-446293302 71^ + 62451523 71^°)] (a; > 12). (72) 



Even for i < 12, the QCM is valid for a number of 
series coefficients, and should coincide with the results 
of the FCM or the "FCM+QCM" . We have used this 
property for checking whether the computer codes are 
correct. 



V. A TECHNIQUE FOR THE SPECIFIC HEAT 

Our strategy for the specific heat in Ref. || was sim- 
ilar to that for the correlation fimctions. That is, the 
most time-consuming part of the FCM is calculated by 
another method instead. If only the FCM is used, the 
series for (Pi^i+i) up to 0[{(3J)^] requires systems with 
e < M/2+1 [M: even) or£ < (M+l)/2 + l (M: odd). In 
Ref. ^, we calculated (Pi^i+i) up to 0[(/3J)^^]. However, 
the FCM was used only for ^ < 11; contributions from 
the required largest cluster, namely, the lowest- and the 
second-lowest-order nonzero contribution from a cluster, 
was directly calculated by a new technique explained be- 
low. 

In fact|n we can use (• • • )o as a moment and define a 
cumulantu without the symmetrization. The generating 
functions of the moments and the cumulants are defined 



as 



you 



)o' (73) 



exchanging positions of uncommutative operators. Note 
that simplification such as Pi,2^\,2 = 1 is not possible in 
[• • • ]o while it is possible in (• • • )o. In expanding the cu- 
mulants using the moments, the ordering of operators 
in (• • • )o must be the same as that in [•■■]o7 for ex- 
ample, [Pi,2P2,3]o = (^,2^2,3)0 - (^'i,2)o(^'2,3)o. This 
cumulant is represented by the cumulant version of an 
amida-diagram. Hereafter, we distinguish between a mo- 
ment amida-diagram and a cumulant amida-diagram by 
bracketing the diagram by (• • • ) and [•••], respectively. 
The free energy is rewritten as 



/3 " /3 ^ m! 

m— 1 



(76) 



(nj'i) (jm,j,„) 

^\P^un■■■P^...,Jo■ (77) 



These are equivalent to Eqs. (|21j) and ( p2D except for 
the subscript "0" replacing "s". Here, however, the cu- 
mulants are generally dependent on the ordering of the 
operators while "s" makes many terms equivalent. 

We can use general property of moments and cumu- 
lants again. If, 



:= lng^J„^'^''='^'^'- ''"''^ (74) 

respectively. Then, (Pi^ • • • Pik,jk)o works as a moment, 
and the corresponding cumulant is define by 



A=0 



Here, the ordering of exchange operators in [• • • ]o is im- 
portant. In other words, the result can be changed by 



\ fc I / 

nexp[AfcP,,,,,]\ {j{e^j,[XiP,,^,,]\ , (78) 
\ fc / \ I / 



, (75) then, the derivative of its logarithm gives. 



= 0. 



(79) 
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When none of ik or jk is equal to ii or ji, this rela- 
tion is obviously satisfied. Furthermore, the cumulant 
is equal to zero also when an amida-diagram is sepa- 
rated into two diagrams by cutting only one point on 
a vertical line as shown in Fig. U^. For example, the 
cumulant versions of the amida-diagrams in Fig. |^(a) 
and (b) are equal to zero. Let us think about the mo- 
ment amida-diagram. Fig. pT|(a), where bonds are lo- 
cated only in the painted regions. Bonds located left 
of ^ are denoted by P^, , and those located right of 
are denoted by Pi^j^- The diagram can be written as 
(Til Piidi Yir Pir,jr)o by exchanging positions of commu- 
tative Pii.ji and Pi^,j^. Let us imagine calculating each 
of YliPiiji and Yl^Pi^j,. by counting cycles, and make 
the product of these two at the end. This final oper- 
ation unites the two cycles relevant to ^ into one cycle 
because Piiji and Hr Pir-J,- sharing only one site 
^. Hence, the number of cycles in the l.h.s. is less than 
that in the r.h.s. "by one" . However, the number of sites 
in the l.h.s. is also less than that in the r.h.s. "by one" . As 
a result, the l.h.s. is equal to the r.h.s. We can make the 



mentioned above, many configurations of bonds make cu- 
mulant amida-diagrams equal to zero. In order to give 
a nonzero contribution, at each vertical line bonds at 
left and right have to appear alternately. We can make 
such diagrams by repeating a enlarging procedure shown 
in Fig. g3| That is, there are two possibilities for next 
bonds. In fact, the value of such a cumulant amida- 
diagram does not depend on those configurations as ex- 
plained in Appendix ^ Hence, we can just multiply 
the value of the cumulants by the probability for such 
nonzero configurations. We have calculated the proba- 
bility by Mathematica using integrals like Eqs. ( ^Sf ) and 



same argument to {Y[i e 



\iPi, 



(a) 



(b) 



FIG. 23: 



The values of the cumulants are shown in Fig. |2j 



Civ n 



(c) example: | ^ | 1 1 = 



-)o. Therefore, Namely, the cumulant is equal to (1 - n~^) when £ is 
an even number, and --(1 — n^^) when £ is an odd num- 
ber. We can prove this by using the mixed expansion 
pivoting on one of Pi, 2 shown in Fig. The first term 
is the moment. The second term is [Pi, 2] times the mo- 
ment of the others. The rest of the terms are obtained by 
separating the diagram at sites from 2 to £ — 1. In fact, 
only the last term survives and every term else cancels 
out one another. When £ is an even number, the sum of 
the first, second and third terms, the sum of the fourth 
and fifth terms, . . . , are equal to zero. When £ is an odd 
number, the sum of the first and second terms, the sum 
of the third and fourth terms, . . . , are equal to zero. The 
details are explained in Appendix 



FIG. 21: Bonds are located only in the painted regions. 
When an amida-diagram can be separated into two diagrams 
by cutting only one point on a vertical line (indicated by an 
arrow), the left and right blocks are independent of each other. 



Here, we must be careful about the periodic boundary 
condition of amida-diagrams. If a diagram satisfies the 
condition above by using the periodic boundary condi- 
tion, the amida-diagram is also equal to zero. We give 
an example in Fig. At a sight, [Pi,2P2,3P2,3Pi,2]o does 
not satisfy the condition above. However, it is equivalent 

to [Pl,2Pl,2P2,3P2,3]0 = 0. 



FIG. 24: 



■Hl(^)-[ffll(ffl)-f^l(H 



FIG. 22: 



Suppose there are two bonds per nearest-neighbor pair 
as shown in Fig. |24. It is a contribution from cluster- 
size £ for 0[{f3jy^~^'^~^] of (Pi,i4_i), or equivalently, 
0[(/?J)^^^~^-'] of the free energy. Because of the property 



FIG. 25: 

We can calculate also the next order by the same tech- 
nique as above. Now, one of the nearest-neighbor pairs 
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FIG. 26: 



APPENDIX A: RELATIONS DERIVED FROM 
THE ISOTROPY 



has three bonds. The value of the cumulant amida- 
diagram does not depend on the position of this triple 
bond as shown in Fig. That is, the cumulant is equal 
to —^(1 — ^) when £ is an even number, and equal to 
— (1 — \) when £ is an odd number. 



VI. SUMMARY 



For example, using r? — \ generators such as 

A'"'^ .= /| (X"'' + X'^") (l<a</3<7i), (Al) 



y^ti ,^ _^^l^'l^X°-(i _ X^'^) (l<a</3<n), (A2) 



We have formulated a new method of high-temperature 
series expansion using the SU(n) Heiscnbcrg model. It 
is designed for efficiently calculating contributions from 
large clusters, which is actually the most time-consuming 
part of a standard method. The net contributions from 
a cluster to high-temperature series start from a certain 
order. We have focused on this property, and considered 
deviation from the lowest-order nonzero contribution. 




FIG. 27: An example of clusters unsuitable for our method. 

Although we have mainly used one-dimensional sys- 
tems to explain the new method, many of the techniques 
can be used also in higher dimensions. In contrast to one- 
dimensional systems, however, the calculation requires 
various clusters, and some of the clusters are unsuitable 
for our method. Figure ^ shows an example. The main 
advantage to use our method is that diagrams can be re- 
duced to smaller diagrams. However, we cannot apply it 
to a cluster in Fig. ^ because it has many sites at which 
more than two bonds meet, and chains connecting those 
sites are too short. Nevertheless, the number of embed- 
dings of such unsuitable clusters in the infinite lattice 
is small, and thus contributions from them are usually 
small. In fact, the number of embeddings of open chains, 
the most preferable clusters for our method, is the largest. 
Hence, the dominant contributions can be calculated by 
our method. 

We expect that similar approaches are possible also for 
other models, namely, for Hamiltonians of lower symme- 
tries. 
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a(pt — 1) 



\1=1 



(2 < a < 7i), 
the exchange operator is rewritten as 

l<Q</3<n 

2<a<n 



(A3) 



,a/3 



1. 



(A4) 



Let us prove that each component contributes equally, 
namely. 



rr — 1 



(A5) 



(i) The Hamiltonian is invariant under re-labeling of 
colors. Therefore, {X"^ X"^) does not depend on 
a or /3. 

(ii) {Xfxf) = {yfyf) for arbitrary a, /3, because 
these X , y operators are x, y components of Pauli 
matrix {x^n/2) regarding only a and (3. 

(iii) We can apply the same argument also for z- 



component, and thus (A;"'';^/'") = {D^' Dy') 
(iv) By a direct calculation using 



for arbitrary a, and 



for arbitrary a ^ f3, we can derive 

(Di'-^D^-^) ^ n[{Xl'Xj') - (A,"Af )) 



for arbitrary a. 
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From (i)-(iv), we obtain Eq. (A5) 
Using Eq. ([A5|), we can derive 



for a ^ /3, z ^ j. 

In addition, for i = j, since the average of |a)(Q;| does 
not depend on a, 



(Xr'^xf) = (|A)(a|) = i^(|7)(7|). (AT) 
Then, using 

n 

V — « i i 

(A8) 



7)(7| = 1, 

7=1 



we obtain 



1 , > 1 



n n 



(A9) 



APPENDIX B: MOMENTS AND CUMULANTS 
FOR CLASSICAL VARIABLES 

In this section, we review fundamental properties of 
moments and cumulants of classical variables Xi {i — 
1,2,--- ). No te that some other properties are written 
in Sec. Ill A. Here, the expectation value with respect 
to the distribution of {xi} is denoted by {•■■)x- The 
generating function of moments is defined as 



That is, the moments are derived by 
{xjXk ■ ■■xi) 



(Bl) 



dXj dXk dXi 



A=0 



(B2) 

where the subscript A = represents that all the A- 
variables are set to zero. On the other hand, cumulants 
are denoted by [■ ■ ■]x and defined as 



\XjXf^ • • ' Xl\x 



d d 



d 



dXj dXk dXi 



ln(e^.-^'^' 



A=0 

(B3) 

Namely, the generating function of the cumulants is writ- 
ten as 



(B4) 



A cumulant is equal to zero if the variables in [• ■ -Jx 
can be partitioned into two groups which are indepen- 
dent of each other in averaging (• ■ • )x. Suppose {xi} is 



independent of {xj}. Let us in advance set Afe to zero if 
Xk is not included in either {xi} or {xj}. Then, 



By taking the logarithm, we obtain 

ln(e^^^'^')x + ln(eS^-^^"^).. 

That is, the generating function of the cumulants does 
not have any term which has a product of Xi and Xj. 
Thus, the corresponding cumulant is equal to zero. 

There are relations between moments and cumulants. 
Moments can be expanded using cumulants, for example, 

{Xi)x = [Xl]x, 

{xiX2)x = [a;ia;2]x + [a;i]x[a;2]2;, 

{xiX2X3)x = [xiX2X3]x + [xiX2]x[X3]x + [xiX3]x[x2]x 
+ [xi]x[x2X3]x + [xi]x[X2]x[x3]x- 

These relations are generally written as 



{Xi ' • • Xl)x — ^ ^ ^ ^ [Xi ' ' ' Xj^x ' ' ' [xk ' ' ' Xni\x- (B5) 



These summations include every partition of x- variables. 
Here, ^ is the number of x-variables in the bracket in 
the l.h.s. Then, the summation denoted by V{£,, C) is 
taken over every partition of ^ elements into ^ groups 
(1 < C < ^). In other words, ( is the number of brackets 
[■ ■ ■]x in the r.h.s. Here, each [■ ■ ■]x includes at least one 

X. 

Conversely, cumulants can be expanded using mo- 
ments, for example, 

[Xi]x = {Xl)x, 
[xiX2]x = {xiX2)x - {Xl)x{x2)x, 
[xiX2X3]x {xiX2X3)x - {xiX2)x{x3)x ~ {xiX3)x{x2)x 

-{xl)x{X2X3)x + '2{xi)x{x2)x{x3)x- 



In general, Eq. ( |B5[ ) is inverted to 

[x.---xe]x = y.y. (-i)^-i(c-i)! 



(B6) 



APPENDIX C: SYMMETRIZATION 

Let us consider expanding W". It includes all the or- 
derings of Pij operators. In other words, it is already 
symmetrized. Therefore, 



Accordingly, 



(CI) 
(C2) 
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Because of a a property of trace, one can cyclically 
rotate variables in the average. In each term of the sym- 
metrization, let us put a certain operator, for example 
Pii^ji^, at the leftmost. Then, k terms give the same or- 
dering. Namely, 



) Jo-(l) ^ *ct(2) J(t(2) ■ ' ■ »o-(fc) ,Jo-(fc) )q 



(2) Jo-(2) ■ ' ■ ^io-(fc-l) J,T(fc-l))o' 



(C3) 



where tr is a permutation of fc — 1 elements. The factor 
k cancels out k in fc! in Eq. (17). That is, the sym- 
metrization of all the operators except one is equivalent 
to the symmetrization of all the operators. Therefore, 
(Pi,jW")o is also already symmetrized, and thus 



(C4) 



APPENDIX D: PROOFS FOR CUMULANTS 

Here, we prove properties of cumulants in Fig. ^. For 
a simplicity, we use the SU(2) notation here. However, 
for general SU(n), one can easily replace spin operators 
below by X"'^, yf, Df ^ in Appendix 

In order to consider Fig. ^(c) and (djT"^ Let us expand 
the l.h.s. of Eq. ( ^6|) using moments, namely, in terms of 
Eq. (B6). Then, let us rewrite the pij operators using 
spin operators in terms of Eq. (0). If a site index i ap- 
pears only once in [•••], every term in the expansion is 
equal to zero because the average at site i is evaluated as 
(sf)o = = {si)o = 0. Consequently, the cumulant 

is equal to zero. 

Note that, in contrast to the Ising model, a moment 
bond-diagram in which an odd number of bonds meet 
at a site can give nonzero value because a product of 
three operators can give nonzero average, for example, 
(sfsj spy sf s|}o. 

Next, we prove Fig. 



IVDl 



IS ex- 



^(e). After Sec. 
plained, this property can be easily proved. Let us think 
about the mixed expansion pivoting on the bond in the 
center. Then, Ci is equal to zero because this mid- 
dle bond always unites two cycles. Furthermore, C2 in 
Eq. (|3^ ) is equal to zero because every cumulant in it is 
equal to zero. 



APPENDIX E: SIMPLE FCM FORMULAE IN 
ONE DIMENSION 



In one dimension, the equations in Sec. Ill E can be 



reduced to simpler forms as shown below. From Eq. (29) 
we can derive 



(Pi-lj-l}^-l 



(El) 



In other words, Eq. (29) is inverted into Eq. (El) us- 
ing full perimeter lattice constants explained in Ref. |l^. 
Combining Eqs. (|3^) and (El), contributions from small 
clusters are canceled out and we obtain a simple formula. 



{Pi.t+xy<( ^'^{Pjj+T^e - {Po-i+^)t-^- (E2) 

j=i i=i 

There are formulae also for the free energy and we have 
used it for calculating the specific heat of the SU(oo) 
limit in Ref. ^. The free energy from the interaction 

part for is := ln(exp(-/37^^''^)}o. The net 



contribution from the £-site cluster is defined as, 

fi--=fe-y2{h + l)fLi,- (E3) 



The total net contribution from clusters smaller than or 
equal to i is defined as 



f'<( E //• 



(E4) 



(=1 



The simplified version of these relations also exists. 
Equation (E3) is rewritten agl3 



fi ~ fe- 2/^_i + ft 



i-2- 



Finally, Eq. ( |E4D is reduced to 

f'<e ^ ft - fe-i- 



(E5) 
(E6) 



Here, f'^^ is correct up to 0[(/3J)2^~i]. Note that /^^ = 

Pint/iVrwhere Pint := -rin(exp(-/37^(p)))o. In^fact, 
Eq. ( p^ ) is reduced to — dPint/dA^ in the ^ 00 
limit. 



APPENDIX F: AN EXAMPLE OF THE TRACE 

Let us consider a five-site system. Then, 

n n 

Tr^""-* 1 = • • ■ {aia20:3a4a5\aia2a3a4a5) = . 



Ql=l Q5=l 



(Fl) 

We show the calculation for P = Pi,4P3,5Pi,2 below. 
Here, P^j- exchanges i-th state and j-th state. 

n n 

= ^ ■■• (Qia2Q3Q4Q5|-P|QlQ2Q3Q4Q5) 

Ql— 1 0:5 — 1 

n n 

= ^ ■■• (aia2Q3Q4Q5|Q4QlQ5Q2Q3) 

Ql— 1 Qs — 1 

n n 



Ql— 1 Q5 — 1 
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Let us consider a condition to give a nonzero contribu- 
tion. Here, ^0:1,04 requires ai — a^, then 5a4,a2 requires 
Oii = ot2, then Sa^^ai requires a2 — ai. These are reduced 
to ai — = a2- In fact, this procedure is equivalent 
to finding a cycUc permutation 1 — > 4 — > 2 ^ 1. By the 
same argument, variables in another cyclic permutation 
3 ^ 5 — > 3 also should be equal to each other, namely, 
Q^3 = Oib- Hence, there are only two independent vari- 
ables in the multiple summation above. Accordingly, we 
obtain 



(5) P - ^2 



Tr^^'P 



(F3) 



In other words, permutation (12345) — > (41523) is com- 
posed of cycles (142) (421), (35) ^ (53), and then, 
the number of independent variables is given by the num- 
ber of cycles, namely, two. 

Let us operate one more exchange operator to P. 
When Pi, 2 is operated, both 1 and 2 belong to a cy- 
cle (142), and this operation breaks this cycle. That is, 
Pi,2P has three cycles (1)(24)(35). In the case of Pi^aP, 
however, 1 belongs to (142), and 3 belongs to (35), and 
these cycles are united, namely, Pi^sP has only one cycle 
(15342). 



APPENDIX G: ANOTHER DEFINITION OF 
QUASI-MOMENTS /CUMUL ANTS 

In fact, a quasi-moment can be defined without a pref- 
actor in Eqs. (|3|) and (p4|). In this case. Fig. |l^(b) is 
equal to a quasi-moment. The generating function for 
this quasi-moments is defined by a derivative as, 



d 



d 



d 



i-PiJ PiA+lPi+l,i+2 ' 



(Gl) 



Consequently, the generating function for the corre- 
sponding quasi-cumulants is written as. 



In^q 



(G2) 



However, we have to integrate the quasi-moment generat- 
ing function j — times to come back to the correlation 
function. Hence, it is difficult to calculate general depen- 
dence on J — I, and thus this definition is not suitable 
for o ur p urpose. This is the reason why we have adopted 
Eq. (IJ). 



APPENDIX H: AMIDA-DIAGRAMS 

First, let us think about moment amida-diagrams. In 
order to give a nonzero contribution, a configuration at 
a nearest-neighbor pair is topologically equivalent to one 
of the two types shown in Fig. Eq. A pair of P^^y can be 
simplified by exchanging subscripts x and y of operators 



between them. We can simplify left and right pairs of two 
bonds to give the r.h.s. of Fig. In the parentheses, we 
schematically represent which sites are included in which 
cycle; sites belonging to a cycle are connected by a line. 



(b) 



i 3 k I 
i 3 k I 



i 3 k I 

-u- 

i j k I 



FIG. 28: A part of amida-diagrams. The pairs of two bonds 
can be simplified. Diagrams in parentheses represents which 
site belongs to which cyclic permutation. 

In the moment versions of amida-diagrams in Fig. |2j, 
we simplify bonds at every other intersite, namely. Pi, 2, 
P3,4, .... The simplified diagrams are made of the pieces 
in Fig. |2^ and special pieces for the ends as shown in 
Fig. When the end pieces are absent, any combi- 
nation of the pieces in Fig. |2^ gives two cycles. When 
£ is an even number, the end pieces are simple vertical 
lines, and do not do anything to cycles. Hence, the dia- 
gram has two cycles and gives n^/n^. When £ is an odd 
number, one of the end pieces is not a simple vertical 
line as shown in Fig. |2^(b), and it unites two cycles into 
one. Therefore, the diagram has only one cycle and gives 
n/n^. The important point is that the amida-diagrams 
of odd £ are equal to those of ^ -I- 1. 



123 4. 5 678 



(a) 



(b) 
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FIG. 29: The simplified diagrams are made of these pieces. 
For a guide to an eye, the pieces are separated by small spaces. 
Diagrams in parentheses represent which site belong to which 
cyclic permutation. 



[N] = (N>-(H)(N)-;^ 



FIG. 30: Explicit calculation of the lowest order. 

Then, we can prove Fig. |2j using a mathematical in- 
duction. In Fig. we explicitly calculates the lowest or- 
der. Next, let us consider the mixed expansion in Fig. |2^ 
for the higher orders. Suppose £ is an odd number. The 
first and second terms cancel out each other by the same 
argument as noted in the context of Eq. (|3^) and Fig. |l^. 
If Fig. U is true for the lower orders, the sum of the third 
and fourth terms, the sum of the fifth and sixth terms. 
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. . . , are equal to zero. As a result, only the last term 
survives. Since the moment amida-diagram in it is equal 
to unity, we obtain —(1 — n^'^). Next, suppose £ is an 
even number. The sum of the first and the second terms 
is 4(1 - A)- The third term is -(1 - ^) times 
Hence, the first, second and third terms cancel out. Fur- 



thermore, if Fig. is true for the lower order, the sum of 
the fourth and fifth terms, the sum of the sixth and sev- 
enth terms, . . . , are equal to zero. Namely, again, only 
the last term survives. Since the moment amida-diagram 
in it is equal to unity, we obtain (1 — n~^). 
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